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By means of time-dependent density-matrix renormalization-group (TDMRG) we are able to 
follow the real-time dynamics of a single impurity embedded in a one-dimensional bath of interacting 
bosons. We focus on the impurity breathing mode, which is found to be well-described by a single 
oscillation frequency and a damping rate. If the impurity is very weakly coupled to the bath, a 
Luttinger-liquid description is valid and the impurity suffers an Abraham-Lorentz radiation-reaction 
friction. For a large portion of the explored parameter space, the TDMRG results fall well beyond 
the Luttinger-liquid paradigm. 
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I. INTRODUCTION 

The dynamics of impurities jiggling in classical and 
quantum liquids has tantalized many generations of 
physicists since the early studies on Brownian motion. 
In particular, the interaction of a quantum system with 
an external environment strongly affects its dynamics . 
Because of this coupling, the motion of a quantum parti- 
cle is characterized by a renormalized mass, decoherence, 
and damping. Polarons 4 , originally studied in the con- 
text of slow-moving electrons in ionic crystals, and impu- 
rities in 3 He 5 are two prototypical examples in which the 
bath is bosonic and fermionic, respectively. These prob- 
lems have been at the center of great interest for many 
decades in condensed matter physics. 

Recent advances in the field of cold atomic gases 6 have 
made it possible to observe and study these phenomena 
from a different perspective and hence to disclose new 
aspects not addressed so far. It is indeed possible to 
accurately tune the coupling between a quantum particle 
and the bath and to modify the many-body nature of the 
bath itself. Furthermore, the dynamics of the dressed 
particle can be studied in real time, thus giving direct 
access to both mass renormalization and damping. This 
problem becomes of particular relevance if the bath is 
one-dimensional (ID). In this case interactions strongly 
affect the excitation spectrum of the bath 7 and therefore 
the effective dynamics of the coupled system. 

The dynamics of impurities in cold atomic gases has 
attracted a great deal of experimental 8 and theoretical 9 
attention in recent years. In particular, Catani et al. 10 
have recently studied experimentally the dynamics of K 
atoms (the "impurities" ) coupled to a bath of Rb atoms 
(the "environment" ) confined in ID "atomic wires" . Mo- 
tivated by Ref. 10 we perform a time-dependent density- 
matrix renormalization group (TDMRG) 11 study of the 
dynamics of the breathing mode in a ID bath of inter- 
acting bosons. 

At a first sight one might think that the problem un- 
der consideration reduces to the study of a particle cou- 
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FIG. 1. (Color online) The width of the impurity breathing 
mode cr(t), in units of fh = (J2/V2) 1 5 (5 is the lattice 
spacing), is plotted as a function of time t (in units of T/2 = 
n/u>2)- Data labeled by empty symbols represent TDMRG 
results corresponding to two opposite values of the impurity- 
bath coupling constant 1112. The strength of interactions in 
the bath has been fixed to ui = 1. The data for U12 = — 0.6 
(squares) are shifted downward by 0.37 £h . The black solid 
lines are fits to the TDMRG data based on Eq. (3). 



pled to a Luttinger liquid 12 . As we will discuss in the 
remainder of this work, it turns out that the noncquilib- 
rium impurity dynamics eludes this type of description. 
This is the reason why we choose to tackle the problem 
with an essentially exact numerical method. Several fea- 
tures of the experimental data in Ref. 10 are also seen 
in our simulations. As we will discuss in the conclusions, 
however, a detailed quantitative account of the data in 
Ref. 10 may require additional ingredients and is outside 
the scope of the present work. Here, we highlight a num- 
ber of distinct signatures of the impact of interactions on 
the breathing motion of an impurity, which are amenable 
to future experimental testing. 
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II. MODEL HAMILTONIAN AND IMPURITY 
BREATHING MODE 

We consider a ID bath of interacting bosons coupled to 
a single impurity confined in a harmonic potential. The 
bath is modeled by a Bose-Hubbard Hamiltonian with 
hopping J i and on-site repulsion U\ > 0: 

Kb = -Ji^SjSi+i+H.c.HE/i^^+^Wii • (1) 

i i i 

Here b\ (k) is a standard bosonic creation (annihilation) 
operator on the z-th site. To avoid spurious effects due 
to quantum confinement along the ID system, we con- 
sider a nearly-homogeneous bath: the external confining 
potential Wi is zero in a large region in the middle of the 
chain (1 < i < L) and raises smoothly at the edges (see 
Sec. A). The local density (fii) — (b\bi) is thus essentially 
constant in a region of length ~ 2L/3. In (almost all) the 
results shown below we fix the number of particles in the 
bath to iVbath = 22 and distribute them over L = 250 
sites, thus keeping the average density to a small value, 
(fii) < 0.1. For this choice of parameters the lattice is 
irrelevant and the model (1) is ideally suited to describe 
a continuum. Indeed, in the low-density limit, Eq. (1) 
reduces to the Lieb-Liniger model 13 (a mapping between 
the coupling constants of the two models is summarized 
in Sec. A). The Hamiltonian describing the impurity dy- 
namics is 

Hi(t) = -J 2 ^(a]a i+1 +H.c.) +V 2 (t)Y / rN i , (2) 

i i 

with Ni — a\di the impurity density operator and i = 
i — «o- Eq. (2) includes a kinetic term and a harmonic 
potential centered at io = (L + l)/2, whose strength V^t) 
depends on time, mimicking the quench performed in the 
experimental study of Catani et al. . Because the on- 
site impurity density (Ni) < 0.15 is low, also in this case 
the model (2) well describes the corresponding contin- 
uum Hamiltonian (Sec. A). In this work we have fixed 
J2/J1 — 2 to take into account the mass imbalance be- 
tween the impurity and bath atoms. For future purposes 
we introduce u)(t) = 2^/V<i{£)J<il%) with u% = uj(t < 0) 
and C02 — u(t > 0), i.e. the harmonic-potential frequen- 
cies before and after the quench, respectively. The full 
time-dependent Hamiltonian H(t) = Hb +Hi(t) +H CO up\ 
contains a further density- density coupling between bath 
and impurity i-L coup i = U12 J2i ^iN. 

The quench in uj(t) excites the impurity breathing 

mode (BM), i.e. a mode in which the width a(i) = 
— 1/2 

[Y^i i 2 (Ni(t))~\ , associated with the impurity density 
(Ni(t)), oscillates in time 14 . This quantity is evaluated 
with the TDMRG. 



III. NUMERICAL RESULTS 

In Fig. 1 we illustrate the time evolution of the im- 
purity width a(t) dictated by %(i) 15 . Different sets of 
data refer to two values of the impurity-bath interaction 
U12 = U12/ J2- Time t is measured in units of T/2, where 
T = 2tt/lu 2 is the period set by the harmonic-confinement 
frequency u>(t) after the quench. The TDMRG results 
(empty symbols) have been obtained by setting u\ — 
U 1 /J 2 = 1 and V 2 {t)/J 2 = 10~ 3 for t < and 10~ 4 for 
t > 16 . 

The black solid lines are fits to the TDMRG data based 
on the following expression: 

^ = i+^£w(n/^-».)- 

i—x,p 

(3) 

where cf> — arccos(Y / l — r 2 /Q 2 ), 9 X = 0, 6 p = tt/2, 
A x = [cr(0)/cr(+oo)] 2 - 1, and A p = [er(+oo)/cr(0)] 2 - 1. 
Eq. (3) is the prediction for the BM width obtained by 
solving a quantum Langevin equation for the impurity 
position operator X(t) in the presence of Ohmic damp- 
ing and a random Gaussian force with colored spectrum 
(see Sec C): 

d?X(t) + 2Td t X(t) + n 2 X(t) = i(t) . (4) 

The three parameters cr(+oo), f2, and T (respectively the 
asymptotic width at long times, the frequency of the 
breathing oscillations, and the friction coefficient) have 
been used to fit the data. The initial width tr(0) is ex- 
tracted from numerical data for the ground-state width 
at t < 0. Note that, in the limit in which the impurity- 
bath interaction is switched off (U12 = 0), a(t) must os- 
cillate at the frequency 2w 2 , since only even states (under 
exchange x o — x) of the harmonic-oscillator potential 
are involved in the time evolution of a symmetric mode 
(like the BM). This implies 57 —> 0J2 in the limit U12 —> 0, 
where Eq. (3) reproduces the exact non-interacting dy- 
namics. 

In Fig. 2 we plot the values of the frequency as a 
function of u 12 and for different values of U\. Several fea- 
tures of the data in this figure are worth highlighting: i) 
the behavior of VI is dramatically different when the sign 
of interactions is switched from attractive to repulsive, 
except at weak coupling, for a tiny region of small U12 
values; ii) the behavior becomes more symmetric with 
respect to the sign of u\2 as the bath is driven deeper 
into the Tonks- Girardeau (TG) regime, i.e., for u\ — > 00 
(in passing, we notice that our results in this limit are 
relevant in the context of the so-called "Fermi polaron" 
problem 17 ); iii) the renormalization of the frequency fi is 
reduced on increasing the strength of repulsive interac- 
tions in the bath U\. It becomes almost independent of 
U12 for u\ > 0.5. 

Fig. 3 illustrates the dependence of the damping rate 
r on 1*12, for different values of u%. Three features of 
the data are remarkable: i) T displays a strong asym- 
metrical behavior with respect to U12 — away from 
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FIG. 2. (Color online) The oscillation frequency of the breath- 
ing mode (in units of 0)2) as a function of the impurity-bath 
Lieb-Liniger parameter U12, for different values of the bath 
dimensionless coupling constant Ui. All the data in this fig- 
ure have been obtained by setting V^ifyjJi = 10 -3 for t < 
and 10 -4 for t > 0. Error bars refer to the fitting procedure. 
Solid lines are guides to the eye. 



the weak-coupling limit; ii) T decreases with increasing 
Mi, saturating to a finite result in the TG limit; and, fi- 
nally, iii) r depends quadratically on u\2 for |ui2| -C 1. 
The non-monotonic behavior of T on the attractive side 
can be explained as following. For U12 ~ the damping 
rate must be small. Increasing |u 12 1 the damping rate 
increases because the coupling of the impurity to the 
bath increases. However, upon further increasing |tti2| 
another effect kicks in. We have indeed discovered (data 
not shown here) that T decreases monotonically with de- 
creasing frequency. As shown in Fig. 2, on the attractive 
side £1 decreases rapidly as |ui2| increases, thereby reduc- 
ing the damping rate. The non-monotonic behavior of T 
does not occur for u\% > because on the repulsive side 
f2 changes slightly with respect to u\2- 

The asymptotic width <r(+oo), shown in Fig. 4, fairly 
agrees with the equilibrium value at the frequency u>2 
that one can calculate numerically. This finding seems 
to suggest that the impurity has nearly "thermalized" 
with the bath over the time scale of our simulations. 



IV. LUTTINGER-LIQUID THEORY AND THE 
ABRAHAM-LORENTZ FRICTION 

We now discuss which features in Figs. 2-4 can (or can- 
not) be explained by employing a low-energy Luttinger- 



FIG. 3. (Color online) Same as in Fig. 2, but for the fric- 
tion coefficient F (in units of 012)- Note that F vanishes 
quadratically for weak impurity-bath couplings ( | ltx2 1 "C 1) 
and saturates to a finite value in the limit ui — ► 00. The 
inset illustrates the dependence of F [in units of u^r] on 
U12 (in units of 1/r). For each value of UJ2 , the tDMRG 
data (green symbols) have been obtained by performing a 
quench corresponding to a value of u>\ = vTO U2- The 
other parameters are: iVbath = 40, L = 600 ({hi) ~ 0.07), 
U12 =0.1, and Ui = 1. The solid line represents the predic- 
tion Fal(<^2 <§C 1/t)/(u)2 t ) — 1/2, based on the Abraham- 
Lorentz model with the value of r corresponding to 111 = 1. 
The dashed line at F/uj^t = 0.44 is the result of a best fit to 
the data. 



liquid description of the bath. 

The Hamiltonian of a single impurity of mass M, de- 
scribed by the pair of conjugate variables (X.P), coupled 
to a bath of harmonic oscillators (the bosonic excitations 
of the Luttinger liquid) with dispersion ujf, = v s \k\, is 1,10 : 

H(t) = ^ + V(X, t) + J2 f^kll% + 912KX) , (5) 

where V(x, t) = Mu> 2 (t)x 2 /2 (the identification with the 
lattice model fixes M = h 2 /(2J 2 5 2 ) with 5 the lattice 
spacing) and gi2 is a coupling constant playing the role 
of U12 in the discrete model. In Eq. (5) 7J (7^) is the 
creation (annihilation) operator for an acoustic-phonon 
mode with wave number k and p(x) is the bath den- 
sity operator. The sound velocity v s is related to the 
Luttinger parameter K of the Lieb-Liniger model by 
Galilean invariance 18 ; K is in turn defined by the rela- 
tion k = K/(ttHv s ), where n is the compressibility of the 
bath 7 . The parameters K and v s , which completely char- 



4 



0.70 


i i i i 


i i i 


0.65 


/ 4' 


* \i 

••*.}} 

. J. 


0.60 


"i 

i 




0.55 


i 






i 


Ul 




i 


• 0.1 

* 0.3 

x- oo 

1 1 1 


0.50 


i 

i i i i 




-0.6 -0.4 -0.2 0.0 


0.2 0.4 0.6 



Ul2 



FIG. 4. (Color online) Same as in Fig. 2, but for the asymp- 
totic value a(+oo) at long times of the width o(t) (in units 
of ^ho). Here the lines are not guides to the eye, but repre- 
sent the equilibrium value <T oq for the width in the harmonic 
potential with frequency u>2 (o~ C q = £h /V^ for un = 0). 



acterize the Luttinger liquid, can be expressed in terms 
of the coupling constants of the model in Eq. (1) (see 
Sec. B). We now observe that the sign of the impurity- 
bath coupling <?i2 can be gauged away from the Hamil- 
tonian (5) by the canonical transformation 7^ — > — 7^. 
This means that if the bath was truly a Luttinger liquid, 
impurity-related observables such as Q and T should not 
depend on 1712 being attractive or repulsive. This low- 
energy description seems to apply only in a tiny region 
around U12 = 0. All the deviations from this prediction 
seen in Figs. 2-4 have to be attributed to physics beyond 
the Luttinger-liquid paradigm. 

The Heisenberg equation of motion induced by the 
Hamiltonian (5) reads (see Sec. B): 



Md 2 X(t) + MujjX(t) + M / dt' t(t,t')d t >X(t') 



-gi2d x p(x,t)\ 



=X(t) 



(6) 



where 



f ( f ' t ') = Ei^ e ^ (t)e ^ (t ' )cos ^( t - t ')] - (7) 



approximate the operator T(t,t') with the following c- 
number T(t - t') = J2k^o c fc fc2 cos [<*>k(t - t')]/{Miul). 
In this limit it is possible to show (see Sec. B) that 
Eq. (6) reduces to a quantum Langevin equation with an 
Abraham- Lorentz (AL) term, which describes the reactive 
effects of the emission of radiation from an oscillator 19 . 
This is a term of the form — MrdfX(t) with a "char- 
acteristic time" r = g 2 2 K j '{it M hvj) . Remarkably, ne- 
glecting the well-known "runaway" solution 19 and keep- 
ing only the damped solutions, we find that the quan- 
tum Langevin equation with the AL term yields an ex- 
pression for a(t) after the quench which is identical to 

Eq. (3) with T = r A L(w 2 ) " 2 = <<a wfr/2. The full func- 
tional dependence of Tal on ui2 is reported in Sec. C. 
Note that Tal is proportional to g\ 2 . This is in agree- 
ment with the TDMRG results shown in Fig. 3 in the 
weak-coupling — > limit. Moreover, Fal is propor- 
tional to K/v* (x K 5 (in our case v s cx K~ x from Galilean 
invariance 18 ) and proportional to The former state- 
ment implies a fast saturation of the friction coefficient 
to a finite value in the TG limit (K = 1). This is in 
agreement with the TDMRG data shown in Fig. 3. The 
quadratic dependence of the damping rate on ll> 2 is also 
well displayed by the TDMRG data, at least for L02T "C I , 
as shown in the inset to Fig. 3. 

It is instructive to compare our findings with the ex- 
perimental data of Ref. fO. The latter show that, in a 
sizable range of interaction strength u 12 , the frequency 
of the breathing mode does not vary appreciably while 
the damping coefficient increases up to T ~ 0.2 w 2 . As 
shown in Figs. 2 and 3, we do observe the same behav- 
ior for ui > 0.2. Moreover, as in the experiment, we do 
see that the width of the breathing mode reduces upon 
increasing Ui2- A detailed quantitative comparison with 
the experiment is, however, not possible at this stage: i) 
one notable difference is that our calculations are car- 
ried out at T = 0, while temperature effects seem to be 
important in Ref. I 0; ii) furthermore, the impurity trap- 
ping frequency in our calculations is considerably larger 
than in the experiment 16 . Extending the current calcu- 
lations to take into account these differences lies beyond 
the scope of this work. 

In summary, we have shown that, in the dynamics of 
an impurity coupled to a ID bosonic bath, a Luttinger- 
liquid description of the bath is applicable only in a very 
small region of parameter space, where the impurity suf- 
fers an AL radiation-reaction friction. Among the most 
striking features we have found, we emphasize the non- 
monotonic behavior of the damping rate and the large 
renormalization of the oscillation frequency for attractive 
impurity-bath interactions. 



withcfc = —gi2[Kv s /(irhL)] 1 / 2 \k\e~\ k ^ 2kc , is the memory 
kernel 1 (k c is an ultraviolet cut-off). 

If the dynamics of the impurity is slow with respect to 
the speed v s of propagation of information in the bath, 
then "retardation effects" can be neglected and we can 



Appendix A: Lattice to continuum mapping 

In our simulations we consider a ID bath of interact- 
ing bosons coupled to a single impurity confined in a har- 
monic potential. The bath is modeled by a Bose-Hubbard 
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FIG. 5. External confining potential for the bath in a lattice 
with L = 250 sites, as defined by Eq. (Al) (we set A = 50 
and W/Ji =0.1). 



Hamiltonian Hb, Eq. (1) in the main text, with hopping 
Ji and on-site repulsion U\ > 0. The external confining 
potential Wi has the following explicit form: 



(W 

2 



1 — tanh ( 



3(i-(A+l)/2) 
2y/(i-l)(A-i) 



ill 

2 

W 



1 + tanh i 



3(i-(2L-A)/2) 



i = l, 

1< % < A, 

A < i < L - A, 

L - A < i < L, 

i = L. 

(Al) 

The parameter A is the length in unit of lattice sites of 
the left and right boundary regions, where the potential 
goes from to a finite value W > 0. The interpolation 
between the two values is as smooth as possible, since the 
potential, as a continuous function of x = iS (5 being the 
lattice spacing), has zero derivatives of all orders at the 
joining points in Eq. (Al). We set A = 50 lattice sites, 
while we used W/J x = 0.1 for L = 250 and N — 22 (see 
Fig. 5), W/Ji = 0.06 for L = 600 and N = 40, where L 
is the length of the (two) lattices used in our simulations, 
and N the number of particles in the bath. In the latter 
case W is smaller since the on-site density rij w 0.07 in 
the middle of the chain is smaller and a weaker potential 
is used. 

With our choice of parameters, the local density (hi) = 
(b\bi) is essentially constant in a region of length ~ 2L/3, 
and kept to a low value (hi) < 0.1 everywhere. At such 
low densities and for not too strong repulsive interaction 
(Ui/Ji < 10), the lattice is irrelevant and the model can 
be mapped to a Lieb-Liniger Hamiltonian 13 describing 
ID bosons of mass m = ft 2 /(2Ji<5 2 ) interacting through 
a contact (repulsive) two-body potential. The upper 
bound for the parameter U\/J\ can be understood as 
follows: as the interaction between bosons is increased, 
the healing length of the Lieb-Liniger gas gets smaller, 
and, when it is comparable with <5, lattice effects be- 
comes relevant. An extensive discussion of this point can 
be found in Refs. 20-22. Other relevant parameters of 
the continuum model are the density n = (hi)/ 8 (where 
(hi) is the on-site density taken in the central region 
where the bath is homogeneous) and the dimensionless 
Lieb-Liniger parameter 71 = mgi/(h 2 n) — Ui / (2J\(hi)) , 



where 171 = U\8 > is the strength of the contact repul- 
sion between bosons. 

The Hamiltonian "Hi (i) describing the impurity is writ- 
ten in Eq. (2) in the main text. Given the low impurity 
density (iVj) < 0.15, also in this case the lattice model 
well describes the continuum Hamiltonian of a particle 
of mass M = h 2 /(2 J2S 2 ) moving in a parabolic potential 
V(x,t) = V 2 (t)(x/8) 2 = Mw 2 (t)x 2 /2 (centered without 
loss of generality at x — 0). The impurity mass is fixed 
to M = m/2 (thus corresponding to a ratio between the 
two lattice hoppings J2/J1 = 2), such to take into ac- 
count the mass imbalance between Rb and K atoms, as 
experimentally done in Ref. 10. 

The total Hamiltonian H(t) = W.B+'Hi(t)+'Hcou P i con- 
tains a further density-density coupling between bath and 
impurity, H CO upi = U\i J2i which in the continuum 

limit corresponds to a <5-function of strength gi2 = U12S. 



Appendix B: Quantum Langevin equation for a 
particle in a Luttinger liquid 



In this section we derive a quantum Langevin equa- 
tion 1,2,23 for a single particle, described by the conju- 
gate variables (X,P), that is coupled to a Luttinger 
liquid by a density-density interaction. The relevant 
Hamiltonian 10 is given by Eq. (5) where the first two 
terms describe the impurity Hamiltonian, the potential 
V(X) = —Muj 2 X 2 /2 representing an harmonic confining 
trap of frequency uj for the impurity, while the third one 
denotes the quadratic Luttinger Hamiltonian, in which 
l\ (7*) 1S t ne creation (annihilation) operator for an 
acoustic-phonon mode with wave vector k and dispersion 
LOk = v s \k\. Apart from the coupling constant 512, the last 
term, which contains an ultraviolet cut-off (fc c ), defines 
the bath density operator p(X) of the Luttinger liquid 7 , 
and is controlled by the parameter K. This, in turn, 
controls the speed of "sound" v s by virtue of Galilean 
invariance: v s = Tmn/ (mK). Repulsive interactions in 
the bath enter the problem through the dependence of v s 
and K on the Lieb-Liniger parameter 71 13 . 

In Eq. (5) the phonon modes couple linearly to the 
particle position, the latter entering through the ex- 
ponential exp (ikX). As a consequence, the harmonic 
excitations of the Luttinger liquid can be integrated 
out leaving an effective dissipative equation (quantum 
Langevin equation 23 ) for the impurity degree of freedom. 
We first switch from annihilation and creation opera- 
tors to (complex) position and momentum, by defining 
x k = y/h/2uj k (% + 7l fc ) and p k = i^tUv k /2(% - j-k), 
with = ihSkk'- Notice that, unlike the usual 

position and momentum, they are complex and their ad- 



joint are x\ 



X-k 



and 



Pk 



-p-k- Using these new 
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variables one gets 
V(X) 



H 



2M 



k>Q 



N 2 



Xk 



C k -ikX 
n e 



+ const. 



(Bl) 



-|fc|/2fcc have 



where the coefficients c k = ~5i2y ^^|fc|e 
been defined according to Ref. 1. Usually, after complet- 
ing the square, some additional terms depending only on 
X are left out and their effect is to renormalize the po- 
tential V(X). Given the peculiar non-linear nature of 
the coupling oc exp(ikX), this does not happen in our 
case: the terms completing the square turn out to be 
independent of X, therefore the potential V(X) is not 
renormalized. 



The Heisenberg equations of motion dictated by 
Eq. (Bl) are given by 

Md 2 X - F{X) = ikc k x k e lkx , (B2) 

fe#0 

8 2 x k + cu 2 x k = c k e~ lkx , (B3) 

where F(x) = —d x V(x). The solution of the second 
equation can be immediately written down: 

x k (t) — x k (0) coswfct + - smuj k t 



+ — dt' sm[u k {t - t')} e- tkx ^ . (B4) 

^k Jo 

Integrating by parts and substituting into Eq. (B2), we 
get 



Md 2 X(t)-F(X(t)) + M / dt't{t,t')d t >X(t') 



k^O 



ikX(t) 



(iCk(O)- 



£k_ -tkX(O) 



P-fc(O) . 

cos uj k t H sin u) k t 

u k 



(B5) 



with a memory kernel function 



f (*,*') = 



elk 2 



^ Mco 2 

fc/O k 



•* k *®e- ik *V ) coau k (t—e). (B6) 



that depends separately on t and t' (and not only on 
t — t'), due to the presence of the two non-commuting 
operators X(t) and X(t') evaluated at different times. 
Upon the substitution x k (0) — c k e~' lkx ^ /oj 2 — > x fe (0), 
the right hand side is just — gi2d x p(x, t)\ x —xit)- ^ n or ~ 
der to see this, simply rewrite x k (0) and p k (0) using the 
corresponding annihilation and creation operators and 
compare to the density-density coupling in the Hamilto- 
nian (5). As discussed in Rcf. 1, one can forget about the 
term oc e lkx ^ in the right hand side of Eq. (B5), since 
it has no effect when taking averages on the equilibrium 
state of the bath coupled to the particle at t = (as a con- 
sequence the stochastic term in the quantum Langevin 
equation (B5) has zero average: (d x p(X(t),t)) = 0). 

When evaluating the noise correlator one runs into dif- 
ficulties since e ^x(t) e -ikx{t') ^ e ik{x{t)-x(t')] _ In the 
limit when the particle has a small velocity with respect 
to v s , we can however perform the following approxima- 
tion: 



JkX(t) -ikX(t') ±iv B \k\(t-t' ) ^ ±iv t \k\{t-t') 



(B7) 



since in this case the product of the first two exponentials 
is assumed to be slowly varying with respect to the third 
one. Therefore in this limit the noise correlator S(i, t') = 



gl 2 ([d x p(X(t),t)?d x p(X(t'),l/)) reads: 
E(t t>) « ^ r -e 



f3huj 



x f coth 
V 2 



cosw(t — t') — ismco(t — t') 



(B8) 



where we passed in the continuum by substituting the 
series over k with an integral over the frequencies u. In 
the same approximation, the memory kernel (B6) enter- 
ing the third term in the left hand side of Eq. (B5) reads 



f(t,f) 



nhvt 



did 

2^ 



e^^Lu 2 cos uj{t-t'). (B9) 



Notice that, when the cut-off to c goes to infinity, this 
kernel basically reduces to the second derivative of a delta 
function: 



t{t,t') u ' c ^" x -Mr8"{t-t'), 



(BIO) 



where r = g 2 2 K / '(ttMHv*) is the proper time scale. 

The relation 5ft = Mhu coth(/3/iw/2)5ft[f (uj)] 

states the fluctuation-dissipation theorem 1,2,23 , and holds 
between the real parts of the Fourier transforms of the 
noise correlator, S(w), and of the memory kernel, T(o;). 
This ensures the consistency of the slow-particle approx- 
imation, see Eq. (B7). Within this approximation, the 
quantum Langevin equation takes the linear form 



Md 2 X{t) - F{X(t)) - MrdfX(t) 



-gi 2 d x p(X{t),t) . 

(Bll) 
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The fact that the position operator appears on the right 
hand side as an argument of the bath density does not 
spoil linearity, since the noise correlator is independent 
of the difference X(t) — X(t') in the slow-particle approx- 
imation. 

Considering X(t) as a classical variable, the fc-sum in 
the third term on the left of Eq. (B5) can be easily per- 
formed, thus obtaining the following classical Langevin 
equation: 



Md 2 t X{t) - F(X(t)) 



9hK 



dt' dtX{t') x 



2-Khvf J 

[S(X(t)-X(t') + ev s (t-t'))] = -g 12 d xP (X(t),t) . 

(B12) 



From here it is apparent that the particle is interacting at 
point X(t) and time t with the phonons (density fluctu- 
ations in the Luttinger liquid) emitted by itself at point 
X(t') in a past time f '. This is similar to radiation damp- 
ing of the motion of a charge particle, where the role of 
the electromagnetic held is now played by the Luttinger 
liquid. The problem of damping due to the emission of 
radiation is a very old one and the classical version of 
Eq. (Bll) has been known in this context for a long time 
as the Abraham-Lorentz equation 19 . 

The quantum version of Eq. (B12) is more compli- 
cated, due to its operatorial character (in this case in- 
deed e *kX(t) e -ikX(t') ^ e ik[x(t)-x(t')]y In the specifiC) 

Eq. (Bll) has been obtained under the assumption (B7), 
thus it is expected to be valid only when the particle does 
not move too fast with respect to the phonons. This can 
be quantified by considering a small value for the ratio v 
between the maximum velocity of the impurity distribu- 
tion when it expands, u)-x$^J o{$$) ( nere ^ho = y/h/Mcu2 is 
the harmonic oscillator length) , and the velocity of sound 
v s . Such ratio is independent of the trap frequency, and 
is fixed only by the initial squeezing through the uncer- 
tainty principle - if the uncertainty in position is cr(0) 
then the uncertainty in velocity is H/(M<r(0)). 

In our simulations this parameter is actually not so 
small. An estimate can be given as follows: the velocity 
of sound v s is a fraction of the Fermi velocity for a gas 
of free fermions at the same density, say v s ~ hnn/m, 
the initial squeezing is er(0)/^h o ~ 0.4 and the density 
niho ~ 1, so 



h _ m/M 
v s a(0) ~ Mv s a(0) ~ ira(0)n 



ho 



1.6. 



(B13) 



So we expect corrections to Eq. (Bll) due to retardation 
effects embodied in Eq. (B5). 



Appendix C: Impurity breathing mode within the 
quantum Langevin equation 

In this section we derive the function used to fit the 
TDMRG data [see Eq. (3) in the main text]. We show 



how it can be obtained starting both from the quantum 
Langevin equation for an ohmic bath 1 ' 2 : 

Md?X(t) + 2MTd t X(t) + M^I(() = |(t) , (CI) 



with noise correlator 



m = mm) 



(C2) 



2Afr 



+00 



dcohuj 



coth cos tut — i sin ujt 



and from the Langevin equation (Bll) in the previous 
section, with a noise correlator given by Eq. (B8) that 
describes a superohmic bath 1 ' 2 . 

The predictions for the oscillation frequency and the 
damping coefficient are different in the two cases. How- 
ever it is possible to treat both equations in the same 
way, because their respective Green functions G(uj) have 
a similar pole structure in the frequency domain. Both 
have two poles at complex frequencies u>± with negative 
imaginary part (r > 0): 



= ± y / n 2 -r 2 -ir. 



(C3) 



These are the physically relevant ones and correspond to 
exponentially decaying solutions. It is essential for our 
purposes that lu + = —uj*_. We point out that Eq. (Bll) 
is of third order and its Green function has a third pole in 
the upper half of the complex plane that corresponds to 
an unphysical diverging solution, also called "run-away" 
solution 19 , which will be discarded in the following. Note 
that, while for the ohmic case f2 is equal to the trap 
frequency W2 , and V is exactly the coefficient that appears 
in Eq. (CI), for the Langevin equation in (Bll) one has a 
more complex dependence on the trap frequency 0J2 and 
on the time scale r. In the specific, the three roots of the 



cubic characteristic equation 



■ or 



w± = ± 



1 



run— away 



3r 



1 f z + z 
3r V 2~ 



[l + z + z- 1 ] , 



read 

. (C4) 
(C5) 



with 



1/3 



_ ±1 = ( 27^ + 2 ± N /(27^ 2 r 2 + 2) 2 - 4 ^ ^ 



Note that u± have always negative imaginary part, while 
^run-away has positive imaginary part and zero real part. 

The symmetry of the physical roots allows us to write 
the solution of both the two Langevin equations as 



X{t) = X(0)e- rt cos v 7 ^ 2 - TH 

V^ 2 - r 2 



(C7) 



dt'Git - t')i{t') . 
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£(t) and G(t) denoting the noise term and the Green 
function respectively. 

Using Eq. (C7), and the fact that = 0, we can 

evaluate some asymptotic averages which turn out to be 
useful to calculate the impurity breathing mode. Some 
of them are null: 

(X(+oo)) = (d t X(+oo)) = ({X,d t X}) t ^ +oo =0, 

while other ones, such as (X 2 (+oo)) and ((dtX (+oo)) 2 ) , 
can be implicitly written in terms of the spectral func- 
tion 1 , whose form depends on the nature of the bath. 
Finally, call E(t) = J* dt'G{t - t')£(t') and note that 



(X(0)E(t)) = (d t X(0)E(t)) =0. 



(C9) 



With all these results at hand, the average of X 2 (t) can 
be calculated at an arbitrary time t. Suppose that the 
particle starts its motion at t = — oo in a squeezed har- 
monic potential. At t = it will have equilibrated with 
the bath, and expectation values are taken on the equilib- 
rium state of the whole system, with the particle at rest 
in the squeezed harmonic potential. Then the frequency 
of the harmonic confinement is suddenly quenched to a 
new value u>2- The width of the impurity density distri- 
bution as a function of time, for t > 0, can be explicitly 
calculated by taking the square of Eq. (C7) on the global 
state: a 2 {t) = (X 2 (t)). Using Eqs. (C8)-(C9), one ends 
up with the following expression: 



a 2 (t) = (E 2 (t)) + - 



-2Yt 



(l 2 (0))cos 2 {^W-TH- 



(CIO) 



with cf) — arcsin(T/f2). It is not necessary to calculate ex- 
plicitly the average (E 2 (t)), since it has to cancel exactly 
the other terms in Eq. (CIO) when (X 2 {0)) = (X 2 (+oo)) 
and ((d t X(0)) 2 ) = ((d t X(+oo)) 2 ) . Therefore the final 
result is 



a 2 {t) =a 2 (+oo)\ 1 + 



-2Tt r 



A x cos 2 (y/W-TH 



+ A p sin 2 y/n 2 - T 2 t 

where the definitions below have been used: 

A = (X 2 (0)) {X 2 (+oc)) 
(l 2 (+oo)) 
A = ((d t X(0)) 2 ) ((cU(+co)) 2 ) 
P fl 2 (X 2 (+oo)) 



(Cll) 

(C12) 
(C13) 



We employ Eq. (Cll) in our fitting procedure as fol- 
lows. A x and A p are fixed to their non-interacting values 



Ho) 



1 and A 



p 



2 (0) 



1, (C14) 



x a 2 (+oo) 

while cr(+oo), Cl, T are used as fitting parameters. The 
initial width <r(0) is provided directly by the numerical 
data. This choice allows reliable fits that smoothly inter- 
polate with the exact non-interacting result: 



a(t) = J cr 2 (0) cos 2 uj 2 t + 



r 4 

1 ho 



4a 2 (0) 



sin u>2t . 



(C15) 



Surprisingly, the same functional form of Eq. (Cll) 
holds for two quite different models. This is due to 
the form of the two relevant complex frequencies [see 
Eq. (C3)] that applies to both of them. However the 
two models provide distinct predictions for the modulus 
O = \lu± I and the damping coefficient T, with some ap- 
preciable differences. For an ohmic bath one has f7 = w 2 , 
so that the frequency is not renormalized in this case. 
Moreover T is a fixed time constant independent of the 
trap frequency. For the superohmic case of Eq. (Bll), 
for LU2T <C 1 it can be shown that Q/lu2 is a function of 
ui2T, the same holds for T/uj2 (see Fig. 6). Notice also 
that uj± always has a non-zero real part, therefore there 
is no overdamped solution for any value of r. 

Finally we emphasize that the functional form given in 
Eq. (Cll) is essential to extract the parameters O, F and 
<t(+oo) from the numerical data. Indeed we have verified 
that the numerical fits performed using the functional 
form 10 



a(t) = cr(+oo) + Aer rt cos(y/n 2 -T 2 t - 0) , (C16) 



with A being the oscillation amplitude (grey lines in 
Fig. 7), are significantly different from the ones that are 
obtained employing Eq. (Cll) (black lines in Fig. 7). For 
example, the fit done with Eq. (C16) often overestimates 
the asymptotic with ct(+oo). 
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FIG. 6. (Color online) Plot of fl/ui2 and T/u}2 as a function of ^/oj^t oc g\2, to allow a comparison with the results in the main 
text. In particular, after expanding Eq. (C4) for u)2T <C 1, we get: Q,/u>2 w 1 — u^ 1 " 2 /2 and T/(c02t) « 1/2. 
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FIG. 7. (Color online) Width of the impurity breathing mode a(t) (in units of £^ = [h/^M^)] 1 ) as a function of time t (in 
units of T/2 = tt/u)2)- Data in different panels corresponds to four values of the impurity-bath Lieb-Liniger parameter u\2, for 
a fixed value of the bath Lieb-Liniger parameter ui = 1. Filled red circles label the tDMRG results. The black solid lines are 
fits to the tDMRG data using Eq. (Cll). The thin grey lines are obtained by employing the fitting function in Eq. (C16). 



